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Abstract We consider a class of pseudodifferential operators, with crossed vector 
valued symbols, defined on the product of two closed manifolds. We study the asymp- 
totic expansion of the counting function of positive selfadjoint operators in this class. 
Using a general Theorem of J. Aramaki, we can determine the first term of the asymp- 
totic expansion of the counting function and, in a special case, we are able to find the 
second term. We give also some examples, emphasizing connections with problems 
of analytic number theory, in particular with Dirichlet divisor function. 
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Introduction 

In Il28l L. Rodino introduced bisingular operators: a class of pseudodifferential oper- 
ators defined on the product of two closed manifolds Mi x M2, related to the multi- 
plicative property of Atiyah-Singer index, see [2|. A simple example of an operator 
in this class is the tensorial product Ai A2, where Ai, A2 are pseudodifferential 
operators on the closed manifolds Mi, M2. Another example, studied in |28|, is the 
vector-tensor product A 1 KIA2. In ||26| . in order to prove an index formula, F. Nicola 
and L. Rodino introduced classical, i.e. polyhomogeneous, bisingular operators and 
defined Wodzicki Residue for this class of operators. The two authors defined the 
residue, via holomorphic families, as in Il9l l25l . For the index of bisingular operators 
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see also the work of V. S. Pilidi l27| and of R. V. Duducava |i5l|6|. In ED, R. Mel- 
rose and F. Rochon introduced pseudodifferential operators of product type, a class 
of operators close to bisingular operators. Bisingular operators are an example of 
operators with vector valued symbols; pseudodifferential operators of this type have 
been meticulously studied, see, for example, Fedosov, Schulze, Tarkhanov [8J and 
the references therein. 

The aim of this paper is to analyze the asymptotic behavior of the counting func- 
tion of selfadjoint elliptic positive bisingular operators. Similarly to the the case of 
SG-calculus [3 1 (see e.g. |7 ,29| for more detail on S'G-calculus), we use techniques 
related to complex powers of operators, i^-function and Tauberian Theorems. This 
strategy, in the setting of closed manifolds, was first used by V. Guillemin 1 14| in 
order to get the so called soft proof of Weyl's formula. 

Here, as in the case of 5G-calculus, it turns out that the i^-function can have 
poles of order two. Thus, using a refinement of Tauberian Theorem due to J. Aramaki 
fl ], the asymptotic behavior of the counting function is determined. The presence of 
a pole of order two of the ^-function implies that the counting functions can have 
asymptotic terms of order A' log A. Such a behavior appears in various setting: mani- 
folds with conical singularities |9 1, 5G-calculus on M" |25 1, 5G-calculus on manifolds 
with cylindrical ends ||20| . See also Gramchev, Pilipovic, Rodino, Wong 110111 II on 
the asymptotic expansion of the counting function in the case of twisted bi-Laplacian. 
Furthermore, in |24|, S. Moroianu studied Weyl's law on manifolds with cusps, with 
an approach similar to the one used in this paper. In a special case, he showed that 
the growth rate of the counting function is A log A. 

We remark that it is not surprising that the ^-function of a selfadjoint elliptic 
positive bisingular operator can have poles of order 2. Indeed, let us consider two 
positive elliptic pseudodifferential operators A,B defined on the closed manifolds 
Mi,M2. From general theory of complex powers of pseudodifferential operators on 
closed manifolds OOl . we know that the i^-function of an operator P of this type is 
holomorphic for fHe(z) < — ^ (« = dim M, m order of P) and it can be extended as a 
meromorphic function to the whole of C with poles of order 1 . As we noticed at the 
beginning, the tensorial product A B is a bisingular operator on My x Mo and it is 
clearly positive and selfadjoint. One can prove the following 



C(A®B,z) = C(A,z)C(B,z). (1) 

If one defines the i^-function using the eigenvalues, equality ([T]) becomes more trans- 
parent. To this end, let {AyjjgN and be the eigenvalues of A and B, respec- 
tively. Then the eigenvalues of A turn out to be {A;7X,};ygpj2. Therefore we have 



C(A,z) = ^A;, $ne(z)<-^; 



C(A0B,z)= X A;m/ = C(A,z)C(B,z), 5He(z) < -max{^, ^}; 
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where ni = dim Mi, n2 = dim M2 and niA^mB are the orders of A and B. Then the 
product structure of i^(A ®B,z) implies that it can have poles of order two. Let us 
now focus on the special case ^ — — zo'- 

C(A,z)=— mc{z)<~zo + e; 

(Z + ZO) 

CiB,z)^-^^+hB{z), mc{z)<-zo + e; (2) 
(z + zo) 

aA^B,z) ^-^^ + ^^f+''fK hAiz)hsiz), mciz) < -zo + e; 

where Ca,Cb are constants that depend just on the principal symbol of A,B, while 
hAjliB are holomorphic functions which depend on the whole symbol of A,B. From 
dU, it is clear that i^(A has a pole of order two. Moreover, we observe that the 

coefficient of the pole of order one depends on the whole symbol of A and B. Finally, 
applying J. Aramaki's Theorem l3.1l from (|2|i one obtains 

zo \ zo Zq J 

(3) 

where 5 > 0. Simple examples of operators A and B for which (O holds are A = 
—Ag + 1, B — ~Agi + 1, where Ag, Agi are the Laplace Beltrami operators associated 
to Riemanniann structures of Mi, M2 respectively. We will extend (O to all posi- 
tive bisingular elliptic operators, expressing the constants in the Weyl asymptotics in 
terms of the crossed vector-valued symbols. 

The paper is organized as follows. In Section[T|we shortly recall basic properties 
of bisingular operators; we refer the reader to r26''281 for more details. Section |2] is 
devoted to the definition of complex powers of suitable bisingular operators; we in- 
troduce the ^-function in this setting and we study its meromorphic extension. The 
main result, concerning the asymptotics of the counting function of self adjoint el- 
liptic positive bisingular operators, is stated in section [3] In section |4] we show the 
connection with Dirichlet divisor problem, which we reconsider from the point of 
view of Spectral Theory. 



1 Bisingular operators 

We start with the definitions of bisingular symbols and bisingular symbols with ho- 
mogeneous principal symbol. In the following, X2, always denotes a bounded open 
domain of K"' . 

Definition 1.1 We define S"''^'"^{Qu£22) as the set of C°°{Qi x i22 x M"i x R"2) 
functions such that, for all multiindex a,-,)?; and for all compact subset K, C X2,, / = 
1,2, there exists a positive constant a2.i3i.fc.Jfi.A^2 

for allx,- e Ki, e E"', / = 1,2. As usual, (^) = (1 + 
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S °°' °°[Qi,Q2) is the set of smoothing symbols. Following ||28| . we introduce the 
subclass of bisingular operators with homogeneous principal symbol. 

Definition 1.2 Let a G gm-mi (X2j ,^2); a has a homogeneous principal symbol if 
i) there exists a„,j^.(xi,X2, (^1,(^2) G 5"i^™2(X2i,i22) such that 

fl(xi,X2,f<^l,<^2) =r'"'a(-Xl,-X2,^l,&), Vxi,X2,^2, V|^i|>l,f>0, 

a — V/i(<^i)fl,ni,- G 5"''"'''"2(X2i,X22), 1^1 cut-off function of the origin. 

Moreover, ai„i, {xi,X2,£,i,D2) G L"j-(i22), so, being a classical symbol on 122, it 
admits an asymptotic expansion w.r.t. the £,2 variable, 
u) there exists a., m2(-^i, -^2,11, fe) e 5'"i''"2(i2i , 122) such that 

a{xi,X2,£,i,t£,2) =t"'^a{xi,X2,£,i.,£,2), Vxi,X2,<^i, V|(^2| > l,f > 0, 
a — ^2[^2)a-.m2 £ 5'"' ''"2 ^'(X2 1,122), 1^2 cut-off function of the origin. 

Moreover, a.^,„^{xi,X2,D\,^2) S (■^1)^ so, being a classical symbol on X2i, it 
admits an asymptotic expansion w.rt. the ^1 variable, 
iii) The symbols a„n^. and a. ^^ have the same leading term, so there exists ami,m2 
such that 



0,^2 - ri(^i)fl™,,m2 e 5""-'''"2(i2i,i22), 

and 

The set of symbols with homogeneous principal symbol is denoted as 5™/ (X2i , i22). 
We will shortly write that the principal symbol of a is {fl„,,..,a. „,2}- 

We can observe a similarity, at least formal, between bisingular symbols with homo- 
geneous principal symbol and SG- classical symbols, see, e.g.. |7,25|. 

We define bisingular operators via their left quantization. A linear operator A : 
C"(X2i xX22) — > C°°(i2i X 122) is a bisingular operator if it can be written in the form 

A(m)(xi,X2) =Op(fl)(xi,X2) 

= fjJn,+n2 I I e"-^'^-^'M^UX2,^U^2H^U^2)d^ld^2. 

If a G5'"i'"'2(X2i,i22) ora G5pr'''"'(^2i,X22), then we write A G L'"'-'"2(r2i,X22) and 
A G Lpr (^1 1 ■^22) respectively. The above definition can be extended to the product 
of closed manifolds; we refer to 1281 for the details of the construction of global 
operators and the corresponding calculus. 
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Definition ll.2l implies that, for every operator A G LjJr' '"^ {£2i , Q2), we can define 
functions a'"' , a'"^ , d"'' ■"'2 such that 

<'(A):ri2i\{0}^C(r22) 

af (A):rX22\{0}^C(^2i) 

(X2,^2) ^ a.,,„2(-«i,-«2,£>l,fe), 
(A) :r*X2i \ {0} X T*Q2 \ {0} ^ C 

,1712 

Moreover, denoting by a{P){x,£,) the principal symbol of a preudodifferential 
operator P on a closed manifold, the following compatibility relation holds 

(T«'(A)(xi,^l))(^2,fe) =CT((Tf (A)(^2,&))(X1,^,) 

= (7'"' (A) (xi ,X2, (^1 , fe) = flmi ,m2 ,-*^2, ^1 , fe)- 

Remark 1 If we consider the product of closed manifolds M\ x M2, then the whole 
symbol is a local object, in general. Nevertheless, similarly to the calculus on closed 
manifolds, it is possible to give an invariant meaning to the functions (HI as functions 
defined on the cotangent bundle, see Il28ll . 

As in the case of the calculus on closed manifolds, it is possible to define adapted 
Sobolev spaces and then to prove some continuity results. 

Deflnition 1.3 Let Mi ,M2 be two closed manifolds. The Sobolev space '"'^ (Mi x 
M2) is defined by 

//'"'■'"2(Mi XM2) - {m e ,r'(Mi XM2) I Op((^i)""(^2)'"^)(«) eL2(Mi XM2)}. 

If MG //"I -'"2 (Ml xM2)then \\u\\,„,.„^, =J|Op((^i)'"i (<^2)'"2)(m) lb- Using the formal- 
ism of tensor product, we can also writqj 

-'"z (Ml x M2) = //"i (Ml ) ® (M2) . 

Similarly to Sobolev spaces H^{M), we have 

i) i/"'i-'"2(Mi XM2) ^ //™i''"2(Mi XM2) is a continuous immersion if m,- > m'-, 
/= 1,2. 

ii) //'"'■™2(Mi X M2) ^ //'"i'™2(Mi X M2) is a compact immersion if m/ > m\, i — 
1,2. 

Proposition 1.1 A pseudodijferential operator A G L'"' ^"2 [m^ x M2) can be extended 
to a continuous operator 

A : //■'■' (Ml X M2) ^ H'-"n '-'"2 (^Mi x M2). 

Furthermore, the norm of the operator can be estimated using the seminorms of the 
symbol. It is also possible to prove the following proposition: 



For definition of see 1321 . 
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Proposition 1.2 Let A e L™' ™2(Mi x M2) be a bisingular operator; if mi <0 (i = 
l,2j, then there exists N G N such that ||A||o,o < sup | £,<7V Pi{<^{x\ 5-^2,^1 ,<^2))|, where 
{pKOIigN f/*^ seminorms of the Frechet space 5'"'''"2(Mi ,M2). 

An operator A G L'"'-'"-{Mi XM2) is elliptic if a,„j,.,fl.„,2,a„,| the three com- 
ponents of its principal symbol, are invertible in their domain of definition. Explicitly: 

Definition 1.4 Let A e L™' "^(Mi x M2); A is elliptic if 

i) c7'«i '«2(A)(vi,V2) ^ for all (vi.vz) G T*Mi \ {0} x rM2\ {0}; 

ii) a"' (A) (vi ) G L'"^^ (M2) is invertible for all vj G r*Mi \ {0}; 

iii) ap{A){v2) G L"J' (Mi) is invertible for all G r*M2 \ {0}; 

where a'"' '"^ (A), a"' (A), CTj" (^) are as in (©. 

In ll28l . it is proved that, if A satisfies Definition 1 1.41 then A is a Fredholm operator 
This property is a corollary of the following theorem: 

Theorem 1.1 Let A G Lj5r '"^ (Mi x M2) be elliptic; then there exists an operator B G 
Lp,."' '"'"^ (Ml x M2) iMc/z that 

AB = ld+Ki, 
BA = U+K2, 

where Id is the identity map and Ki , A'2 are compact operators. Moreover, the symbol 
ofB is b = {(J™' iA)-\ap{A)-^}. 

The proof of Theorem ll.ll is an easy consequence of the global version of the follow- 
ing lemma: 

Lemma 1.1 Let A G L'"i-'"2(X2i x Q2) andBe L'"'p'"'2 (X2i x Q2), then 

{(«°^)mi+m'i,-'(«°^).m2+"4} = H'l ■• ,■ ' «-:™2 % ^-.m^} 

w/iere 

(flo^j fo)(xi,X2,Z)i,<^2)(M) = fl(xi,X2,Z)i,(^2)ofe(xi,ji:2,Di,^2)(M) Vm G C"(i2i), 
{ao^^b){xy,X2,t,i,D2){v) ^a{xi,X2,B,i,D2)ob{xi,X2,t,i,D2){v) Vv G C~(i22). 

In first row the composition is in the space L°°{Qi) of pseudodifferential operators 
on Q.\, in second row, it is in the space L"'{Q,2). 

2 Complex powers of bisingular operators 

In this section we define complex powers of a subclass of elliptic bisingular operators. 
The first step is to give a suitable definition A -elliptic operators w.rt. a sector of the 
complex plane A . 

Definition 2.1 Let A be a sector of C; we say that a G 5™/''"' (Mi ,M2) is A-elUptic 
w.rt. A if there exists a positive constant R such that 
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i) 



((!'"' ■'«HA)(yi,V2) -A) ' e5-'«''-'«2(Mi,M2), 
for all |v,| >R,i= 1,2, and for all A e A. 

ar'(A)(vi)-AldM,eC(M2), 
is invertible for all | vi | > /? and for all A G A . 

af (A)(v2)-AldMi ^L^'iMr), 
is invertible for all | V2 1 > /? and for all A G A . 

In the following, in order to define the complex power of A, we assume that A is a 
sector of the complex plane with vertex at the origin, that is 

A = {z e C I arg(z) G [K-e,-K + d]}. 



ii) 



iii) 



^.^arg = n — 6 




— arg = —K+d 





Lemma 2.1 Leta€ 5""i''"2(i2i,i22) be A-elliptic. For all Ki C 12,-, = 1,2, there exist 
Co > 1 and a set 

i2^,,^, := {z G C\ A I l(^i)'»i(^2)'"^ < \z\ < co(^i>"" {^iD (6) 

^"^^ Co 

spec(a(xi,X2,^i,^2)) = {A G C I a(xi,X2,^i,fe) - A = 0} C X2|j_^j, 

VjCiGA,^,-GK"'; 

|(A-a.„«,(xi,X2,§i,^2))"'|<C(|A| + (^ir(^2r)-i, 

|(a,,,-Ald^J-^|<C(|A| + (^ir(^2r)-i, 

|(a,,,-Ald^,)-i|<C(|A| + (^ir(^2r)-i, 
Vx,- G G M"',A G C\X2^j,5^,j = 1,2, 

w/jcre (ami,- ~ ^'^^i ) ' stands for the symbol of the operator {ami,- i^^ i^2i^\ 1^2) 
A Idi3 an<i similarly for (a.,m2 ~ Id^j^ ) ^- 
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The proof of Lemma lZTl is essentially the same of the one of Lemma 3.5 in ll22l . 

Next, we prove that, if A A -elliptic, then we can define a parametrix of (A — 
A Id). Actually, we prove that, for |A | large enough, the resolvent (A — A Id)^^ exists. 
Restricting ourselves to differential operators, we could follow formally the idea of 
Shubin (|31 1, ch. II) of parameter depending operators. For general pseudodifferential 
operators, it is well know that this idea does not work, see [! 12|. 

Theorem 2.1 Let A e Lpr ''"^ {Mi x M2) be A-elliptic. Then there exists R e R+, such 
that the resolvent (A — Ald)^' exists for A G Ar = {A G A | |A| > R}. Moreover, 

||(A-Ald)-'||=0(|Ar'), AeA«. 

Proof First, we look for an inverse of (A — A Id) modulo compact operators, that is 
an operator fi(A) such that: 

(A-A)oB(A) = Id+7?i(A), A7;i(A) e L"'-"'(Mi XM2), 
B(A)o(A-A) = Id+/?2(A), A/?2(A) eL"''"'(Mi XM2), ^'^^ 

uniformly w.r.t. A G A. In order to find such an operator, we make the principal 
symbol explicit: 

fl-A =psym(fl)-A+c, c e 5'"'"'''"2-i(Mi,M2), 

where psym(fl) = y/ia,,^,. + V^2a-,m2 ^ WiV2<^m,.m2- As we have noticed in Theorem 
11.11 we can write the symbol of the inverse (modulo compact operators) of an el- 
liptic operator In this case we need to be more careful because of the parameter A. 
Following the same construction as in Theorem ll.il we obtain 

/7(A) = { («" (A) - A Mm,)-' , (cjf (A) - A Mm.)-' }. (8) 

The above definition ^ is consistent in view of the A -ellipticity and of the following 
relation 

C^(((^r(^) -■^ Mm2)-'(xi,^i))(x2,.^2) = {a„^^,„2-Xy\xl,X2,^l,^2), 
a(((72"^(A)-A IdMi)-'(x2,fe))(jci,^i) = (flmi,,„2 -A)-'(xl,x2,<^l,<^2)• 
Using the rules of the calculus and Lemma ITTl we can check that B(A) satisfies 
conditions (|7]). By parameter ellipticity, we get that ^i(A) and R2{^) are compact 
operators for A G A, namely 

(A-Ald)oB(A)=Id+/?i(A), 
(A-Ald)oB(A) =Id+/?2(A), 

A/?i(A),A/?2(A) e 5"^-"^ (Ml XM2) uniformly w.r.t. A G A. Sofi(A) is a parametrix 
and its symbol b{X) has the following form 

^(A) = -(flmi,m2(xi,X2,^i,^2)-A)-Vl(&)r2(<^l) 

+ (fl„|,. - A IdM,)-'(jci,JC2,^i,fe)ri(^i) 

+ (a.,m2 -■^ WMi)-'(jCl,JC2,^l,fe)'/2(&), 
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where {om^,- — A 1(1^2) H-'^ii-^Zi^iifc) is the value of the symbol of the operator 
(fl„,i^.(xi,X2,^i,D2)-A IdM2)"' at (JC2,&), and similarly for (a.^,,,, - A Mmi)"'. Fur- 
thermore, denoting by ri (A) the symbol of /?i (A), we easily obtain 

ri(A) = (fl-psym(fl))o/7(A) + (psym(fl)o/7(A))-l, (10) 

hence f"i (A) e 5'^''^' (Mi ,M2) is the asymptotic sum of terms of the type 

dl^gD^lDT-b{X) g e (Ml ,M2). 

Clearly (flmi.m2(-^i ;-*^2, , &) ~ A)^' = (9(|A|^'). By the theory of pseudodifferential 
operators on closed manifolds, the same property holds for the symbols of the opera- 
tors (fl„,,,.(xi,X2,^i,D2) - A IdM2)"^ and {a.^m^{xi,X2,Di,£,2) - X Umi)'^ and their 
derivatives. Thus ri (A) = (9(|A|^^), as a consequence of the calculus. By Proposition 
11.21 this implies ||/?i||/^2 = (9(|A|^^), and the same is true for the operator /?2. So we 
can choose A large enough such that Ri,R2 have norm less than 1 . In this way, using 
Neumann series, we prove that (A — A Id) is one to one and onto, therefore invertible, 
by the Open Map Theorem. Again, by Neumann series, we obtain B{X) such that 
(|9]l is fulfilled with Ri,R2 smoothing and still with norm (9(A^'). Now notice that 
A [B(A) -B(A)] e 5-'"i-i.-'"2-i for all A G A. Furthermore, if we multiply both 
equations in ©by (A — Aid) ^' we obtain 

(A-Ald)-i =B(A)+B(A)7;i(A)+7;2(A)(A-A)-i/;i(A). 

Hence ||(A-Ald)"'l| =C>(|A|"') andA2[(A-A)"'-B(A)] is a smoothing operator 
in L~°°~°°{Mi X M2), uniformly w.rt. A. 

In order to define complex powers of an elHptic bisingular operator, we introduce 
some natural assumptions. 

Assumptions 1 7. A e 5'"i''"2(Mi,M2) is A-elliptic. 

2. o{A)r\A—^(in particular A is invertible). 

3. A has homogeneous principal symbols. 

Remark 2 If we consider a A-elliptic operator A G Lpr' ''""(Mi x M2) with m,- > 
(/ = 1,2), then C7(A) is either discrete or the whole of C, because the resolvent is 
a compact operator (fSTl, Ch. I). Since by Theorem 12.11 we know that for large A 
the resolvent is well defined, it turns out that the spectrum (t(A) is discrete. Then, 
modulo a shift of the operator, we can find a suitable sector such that Assumptions [T] 
is fulfilled. 

Definition 2.2 Let A be an operator fulfilling Assumptions[T] Then, we can define 

A,:=— / l'(A-Xld)-^dX, me(z)<0, (11) 

271 JdA + 

where Ag = A U {z G C | |z| < e}. 

The Dunford integral in (fTTI ) is convergent because ||(A — Ald)^'|| = (9(|A|^') for A 
large enough. As usual, we next define 

A-':=A,_^oA^ mc{z-k)<0. 
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Remarks In Assumptions [T] we require A n C7(A) — 0, therefore in particular the 
operator must be invertible. It is possible to define complex powers of non invertible 
operator as well, provided the origin is an isolated point of the spectrum, see, e.g., 
IHJ. For example, one can define the complex powers of A = —A ® —A on the torus 
X §\ even if A has an infinite dimensional kernel. 

Theorem 2.2 If A G (Mi ,M2) satisfies Assumptions\l\ then A- e L^'-^z-YMi x 

M2) and it has homogeneous principal symbol. Moreover, by Cauchy TheoreniQ 

Y, 

(12) 



a^ 


— {^mi ,m2 






^■,m2Z 





Proof As a consequence of a general version of Fubini's Theorem, denoting by a^ 
the symbol of A^, we obtain 



X\a~XU)-')dX, $Ke(z)<0. 

2k Jd+Ae 

where (a — Ald)^' is the symbol of the operator (A — Ald)^'. By Theorem 12.1 



we know that 
symbols, we have 



(A - Aid) ' -B(A) e L °°- °°(Mi x M2) so, up to smoothing 



(13) 

r~{b{X))dX, 



I 

2n J Of 



where ^2,^,,,^, is as in Lemma IZTI and the second equality in ( fTST l follows by Cauchy 
integral formula. Now, by Lemma IZTI and by the explicit form of ^(A), we get A-' G 
L"'i^''"2^(Mi X M2). In order to show that A' has homogeneous principal symbol, we 
write 

(^(A)) =v/i((7™' (A) - AldM^)-' + V/2(CT"'2(A) - AldM.)-' 
-V/ir2(ff'"'-'"2(A)-A)-'+c(A), 

where Ac(A) G S-"'^-^--'"^-^{Mi ,M2), VA G A. We split integral in O so that 

2kWA->i((T'«'(A)-AldM,)-' (14) 
+ 2^ Id^A, ^'Wiicy'"' (A) - A Mm, r'dX (15) 
- 2^ Sd^A, Vi V/2 (CT"' '"^ (A) - A ) - IJA (16) 

+ 2k.W^^^c(^)«'^- (17) 

The theorem follows from theory of complex powers on closed manifolds for the 
integrals (fT4l i and (flSl l. and from Cauchy Theorem for integral ( fTSI l. Finally, we notice 
that integral ( fTTl i gives a symbol of order {niiz — 1,W2Z — 1). 



^ In equation )12t a~i^. .,a~„j^-^,a^^-.jjj^^ represent respectively (j"''"'(A~),c7™^''(A^),(T'"'~'"'2'(A~), while 
{a„^.y, (a. jn^)' are complex powers of the operators CT™' (A), CT™^ (A) and (a,„^,„^y is the complex power 
of the function fj'"' '"2 (A) . 
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We now introduce the function CQ^iZ) of an elliptic operator that satisfies Assump- 
tions [T] The proof of the following property is similar to the case of compact mani- 
folds (see|31|,ch. II). 

Proposition 2.1 LetAe L'" ' '"^ (Mi XM2), nii > 0, i = \ ,2, be a selfadjoint operator 
satisfying Assumptions^] Then we have 

A\u) = £!;(/;■,«), 

where {AjjygN the spectrum of A, and {./}};gn ^'re the corresponding orthonormal 
eigenfunctions. We define 

C(A,z):=^A;, jne(z)<min{-^,-^}. 
The definition of C,{A,z) in the general case is the following: 

Definition 2.3 Let A G L"'i "'2(Mi XM2) be an operator satisfying Assumptions [T] 
then 

C,{A,z) :— / KAz{x\,X2^x\,X2)dx\dx2^ 9le(z)mi < — ni,9le(z)m2 < — n2, 

where Kaz is the kernel of A-. The integral is well defined if $He(z)mi < — ni and 
£He(z)m2 < — ^2 since, in this case. A' is trace class. 

Theorem 2.3 KAz{x\,X2,yi,y2) is a smooth function outside the diagonal. Further- 
more, Ka~{x\,X2,x\,X2) restricted to the diagonal can be extended as a meromorphic 
function on the half plane {z G C 1 5He(z) < min{ — + e} with, at most, poles 

at the point Zpole = ™ri{^^7 ^^}- ^^^^ po/e can be of order two if ^ — other- 
wise it is a simple pole. 

Proof By definition, the kernel of A' has the form 

Kaz{xi,X2,xuX2) = [ [ a'{xuX2,^u^2)d^id^2- (18) 
First, let us consider the case ^ > Then, if D\t(z) < A~ G L"'i' "'2'(Mi x 

' mi m2 mi ' ^ ^ 

M2) C L^"i^'^'^"2-£(Mi X M2); hence it is trace class and the integral of the kernel is 
finite. We can write a' — af„^^ . + aj, aj G 5'"''^'^™2Z(Mi ,M2) and we have then 

Kaz{x,x) ^ ^ [ f {al, +a;)d^id^2 

+ 7;:: — ^ — — / / (flm,- +al)dtid^2- 

The second integral in ST9[ is an holomorphic function for *He(z) < — ^ + £ since 
we integrate w.rt. the ^1 variable on a compact set. The same conclusion holds for 
the integral of af. on the set { ( ^ 1 , ^2 ) | | ^ 1 1 > 1 , G R'" } because it has order (m 1 z — 
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l,m2z). In order to analyze the integral of a^,^ ., we switch to polar coordinates and 
we obtain 

■'^^I'^fc = / / a„, ^,.deid^2- (20) 

|^,|>l ' " OTlZ + ni JR"2 Js"i-i 

Clearly ( |20] | can be extended as a meromorphic function on {z e C | 9\c{z) < + 
e}, and, moreover, by ( fTSl i. we get 

lim (z+—]Kaz{xi,X2) ^ - /o^x„\„, — / / a,„"[} deid^2- 



ILL ^ 111 
principal symbol. First we write 



The case ^ < ^ is equivalent, by exchanging mi and OT2. 

The case — = — is a bit more delicate, since we have to analyze the whole 



The definition of principal symbol implies that the second term in (1211 1 belongs to 
S'"^--^-""-~-\MuM2), hence the second integral is well defined for fHe(z) < + £ 
and holomorphic for 9le(z) < — ^ + £■ Now we have to analyze the integral of the 
principal symbol. Splitting M"' x ife"^ into the following four regions 

{(^1,^2) I i^ii < T,i<^2i < n, {(fe,^2) I i<^ii < > n, 
{(<^i,&) I i^ii > T, ifel < n, I i<^ii > T,i^2i > n, 



one gets 



"1 7R"2 
(mi+m2)2+«i+n2 



2 ■ 



imiz + ni){m2Z + n2) Js-'i-i Js«2-i "'«iz.'«2z^^l^^2 



^m2Z+n2 



m2Z + n2) i|^i|<T JS'"2-i ' 



(22) 



■miz+ni 



{miz + ni) J\^. 



(«V)i2Z ~ «miz,m2z)'^^l'^^l 



(OT2Z + n2) i|^i|>TJS"2-l 

where h{z) is an holomorphic function for [He(z) < Zpoie + £• The evaluation of the 
integrals in (l22l) are similar to Proposition 3.3 in li25J . and Theorem 2.2 in |i3J. This 
concludes the proof. 
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Since Mi ,M2 are closed manifolds, Theorem l2.3l implies the following: 

Corollary 2.1 Let A e L'"^-'"^{Mi x M2) be an operator satisfying Assumptions\J} 
then tl^{A^z) is holomorphic for 5He(z) < min{ — — ^} and can be extended as a 



ni n2 
mi ' m2 

Laurent coefficients of^(A,z) at z = Zpoie — '"'^ 



meromorphic function on the half plane 9^e(z) < min{ — — ^} + £■ Moreover, the 



_ "1 



lim ( z + — ) C(A,z) = -i:^^^;^;^ [[ II, <■ ddid^2, (23) 
"■ ' nil J {27tj"i+"2mi JJM1XM2. ' 



if^>^ 



li'q, (^ + l^)C(A,z) = - ,,,,„,,„,„., //., L. L_,«J?«'M<^i, (24) 



m2 mi 



rei^(A) = lim {z + iyi;iA,z) = 

lim(z + l){aA,z)-'-^^)=-TRi,2{A) + TRe{A), (26) 



w/iere 

TRi.2iA) :-- 



- — ^— — lim(— // / / (a,„, .y' ~res^(A)logT) 
/ ^ lim (— // / / (a. -rei2(A) log t) 



anc/ 



r/;e(A) := ,^ r/ / ,/ aj m ^oga,ni,>n2dQ\dd2, (28) 

(27r)"l+"2(miTO2) iM, XM2 JS"i-' Js«2-i ^ 

m] m2 

In dZTl l. {ami, y and {a.^m^Y are the symbols of the complex powers of the operators 
cimi.Xxi ,X2,^i ,£>2) and a. „,^{xi ,X2,Di ,^2) ■ In order to obtain the terms in (l26b . (l27b . 
(l28T l. we notice that the constant T in (l22l l is arbitrary and the Laurent coefficients 
clearly do not change if we change the partition of M"' x R"2, therefore we can let T 
tend to infinity. In this way both the fourth and fifth integral in ( |22] | vanish, due to the 
continuity of the integral w.rt. the domain of integration. The evaluation is similar to 
the proof of Theorem 2.9 in 13] and of Proposition 3.3 in ||251 . 
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3 Weyl's formula for bisingular operators 

In this section we study Weyl's formula for positive selfadjoint bisingular operators 
that satisfy Assumptions[T] In the sequel we use the following Theorem, proved by J. 
Aramaki 1 1 1 : 

Theorem 3.1 Let P be a positive selfadjoint operator satisfying Assumptions [7] If 
l^{P,z) has the first left pole at the point —zo fl«fi0 

j~i 0-1)' V*/ z+zo 

extends to an holomorphic function on the half plane {z G C | 5Ke(z) < — zo + e}, 
then, setting 

NFiX)= 1' 

we have 

*(A)~f^(^)""(i^)|„., + 0(A-.-), 

for a certain 5 > 0. 

Theorem 3.2 Let A £L"'i'"'2(Mi XM2) be a positive selfadjoint bisingular satisfying 



AssumptionsU] then 

( C,X'logiX) + C[X' + OiX'-'^)for ^ = ^ = / 
Na{X)^Ic2X^-+0{X^2-'^) forg>^ , A^-, (29) 

for certain 5; > 0, i ~ 1,2,3. The constants C\ , C[ , C2 , C3 depend only on the principal 
symbol of A. 

Proof We use J. Aramaki's Theorem l3.2l which gives the asymptotic of Na (A ) know- 
ing the first left pole of the zeta function. As a simple application we get ( |29] l with 

^^=(2;r)'n+i(„i,„2)/lxMjs".-'L-.("'"""'^)"'^^'^^^ 



" (2;r)«i+«2 («2 »Ji ) JJm, xM2 is"i M-^ ' ' 

r«,,w + r«.w __^,. , J ,K„,rVM«*; (30) 

/ WiWo IJ.^H/I- /sn|-l /sni-l 



Z nin2 jJMi xM2 JS"i-' JS"2 



1 JS"2- 



(27r)"i+«2„2 JJM1XM2JWL" 



-1 



52 • 



^ The Aramaki's Theorem actually requires another assumption on the decay of F(z)f on vertical 
strips. In this case such condition is fulfilled, in view of the relationship between ij-function, heat trace and 
gamma function, see I13I21I . 
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Remark 4 In this paper we are focused just on bisingular operators with homoge- 
neous principal symbol, since our aim is the study of the corresponding Weyl's for- 
mulae. We do not introduce classical bisingular operators and we do not investigate 
the relationship between the poles of the ^-function and Wodzicki Residue defined in 
126 1 . Nevertheless, extending the results of section|2]to classical bisingular operators, 
one can prove that, for a classical elliptic bisingular operator A € L"'''"'^{Mi x M2) 
that admits complex powers, 

Wres(A) := mim2 lim(z - l)^i^(A,z), 

z^l 

where Wres(A) is the bisingular Wodzicki residue defined by Nicola and Rodino in 

He). 



4 Examples 

First we consider the operator A — —A ® —A on the torus §' x We clearly have 
a{A) — {n^rrp-}i„ ,„)eN2. Hence the spectrum is countable and consists only of eigen- 
values. The eigenvalue {0} has an infinite dimensional eigenspace, while all other 
eigenspaces have dimension four Therefore we get 

A^a(A)= 4. (31) 

Let us define the function : N ^ N so that d{h) is equal to the number of ways 
we can write h = m-n, with m, n natural positive numbers or, equivalently, it is equal 
to the number of divisors of h. This function is often called Dirichlet divisor function. 
By a simple computation, we obtain 

A^a(A2)=4D(A) =4 (32) 

Noticing that tl,{A) — Atl,R{2z)C,R{2z), where C,r{z) is Riemann zeta-function, we can 
easily find the coefficients of the asymptotic expansion and we have 

£)(A) -Alog(A) + (27-l)A+0(A>-^), A ^00, (33) 

where 



7:= lim 



(34) 



is the well known Euler-Mascheroni constant. The asymptotic expansion ( l33T l is well 
known (see [17] for an overview on Dirichlet divisor problem; see also [18 , 19|). It is 
still an open question to understand the behavior of remainder. In ifTSl . G. H. Hardy 
proved that 0(A ? ) is a lower bound for the third term. The best approximation, found 
by M. Huxley in 1 16|, is C»(A'(logA)''), where 

131 18627 

c:= 0,3149038462 d := + 1 - 3,238822115. 

416 8320 
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The conjecture is that the remainder is (9(A ' ). 

It is nevertheless interesting to investigate the Unk between Dirichlet divisor func- 
tion and the above resuhs on the spectral properties of a suitable operators. Let us 
notice that in dSTT i we have a slight abuse of notation, since N{k) was only defined 
for positive operators. In this case A = —A (g) —A is non-negative, but has a non trivial 
kernel. In other words we actually consider 

where i\erA is the projection on the kernel of A. This definition is compatible with 
the definition of complex powers of non invertible operators in |4|. The variant of 
our theory to such a setting, which is possible, will be not detailed here. Rather, let 
us now consider the operator := {—A +c) (g) {—A +c), c > 0, defined on the torus 
§' X Clearly, Ac satisfies Assumptions [T] thus we can apply Theorem 13.21 It is 
easy to see that the eigenvalues of Ac are {(n^ +c){m^ +'^)}(n m)GN2' ^ach one with 
multiplicity four. Hence 

N{Ac;X^) = 4 tl{ real numbers of the form (n^ + c) (m^ + c) \ 
{n^ + c){m^ + c) < X, n,m e N}^4Dc{X). 



By Theorem IZ21 we know that a"''"^(Af ^) = {a^-^{Ac)y^ so the constant Ci in 
(|30] | can be easily evaluated 



Ci = -- 



1 1 



r(2;r)2 4 = 2. 



2(2;r)2' 

Since in this case we know the eigenvalue of the operator, TR{Ac) turns into 

1 



(35) 



TRiaiAc) 



2 lim 



4 lim 

T->oo 



,-=-[T] (C + /2)I 



21ogT 



J=0 (c + /2) 



--logT 



(36) 



= 47,.. 



We have named this constant jc because of the link with the usual constant of Euler- 
Mascheroni 7 in ( l34l i. Notice that, letting c tend to 0, Jc goes to +°°; while, if c tends 
to infinity, Yc goes to —00. Finally, we obtain 



D,(A) = -A?(A,;A') 



'Alog(A) + (27,-l)A + 0(Ai-^), A 



(37) 



In this case, knowing exactly the eigenvalues of the operator, we can check our es- 
timate with a numerical experiment. We have checked dJTb for Dc{X) with A = 
10.000.000. In the second column of the Table [T] there is the estimate of the coef- 
ficient of first term of the asymptotic expansion obtained with the software Maple 15, 
in the third the coefficient obtained by ( l37b . and in the fourth the error. We can notice 
that the error increases with c. This is not surprising, since ( |35] | does not depend on c. 
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Table 1 1 st. term approximation 



c 


1 St. term with Maple 


1st. term in (371 


error 




2 


1,024846785 




0,024846785 




3 


0,9916281891 




0,008371811 




4 


0,968979304 




0,031020696 




5 


0,951859819 




0,048140181 




6 


0,938130598 




0,061869402 




7 


0,926687949 




0,073312051 




8 


0,916888721 




0,083111279 




9 


0,908326599 


1 


0,091673401 




10 


0,900728511 




0,099271489 




11 


0,893902326 




0,106097674 




12 


0,887707593 




0,112292407 




13 


0,882038865 




0,117961135 




14 


0,876815128 




0,123184872 




15 


0,871972341 




0,128027659 




16 


0,867459966 




0,132540034 




17 


0,863235614 




0,136764386 




18 


0,859265437 




0,140734563 




19 


0,855520776 




0,144479224 




20 


0,851977951 






0,148022049 




Table 2 2nd. term approximation 








c 


2nd. term with Maple 


2nd. term in (371 




error 




2 


0,40048285 


0,401484386 


0,001001536 


3 


-0,13493765 


-0, 1339381238 


0,000999526 


4 


-0,499994550 


-0,498993281 


0,001001269 


5 


-0,775928050 


-0,774926584 


0,001001466 


6 


-0,997216950 


-0,996213733 


0,001003217 


7 


-1,181650650 


-1,180647904 


0,001002746 


8 


-1,339595550 


-1,3385899520 


0,001005598 


9 


-1,477600650 


-1,476592538 


0,001008112 


10 


-1,600067350 


-1,599058126 


0,001009224 


11 


-1,710092450 


-1,7090842470 


0,001008203 


12 


-1,809939750 


-1,808931287 


0,001008463 


13 


-1,901308850 


-1,9002985710 


0,001010279 


14 


-1,985505550 


-1,9844949070 


0,001010643 


15 


-2,063562050 


-2,0625496430 


0,001012407 


16 


-2,136292950 


-2,1352865400 


0,001006410 


17 


-2,204381450 


-2,2033750580 


0,001006392 


18 


-2,268373150 


-2,2673662890 


0,001006861 


19 


-2,328729950 


-2,3277195600 


0,001010390 


20 


-2,385833550 


-2,3848212840 


0,001012266 



In order to make the error smaller, we should increase the number of digits at which 
we truncate the series Dc{X). In Table |2] we analyze the coefficient of the second 
term. In this case the error is essentially independent of c, this is due to the fact that 
( |36] | does depend on c. 

Our spectral approach to Dirichlet Divisor function suggests that others Weyl's 
formula techniques (e. g. Fourier Integral Operator) could be useful to attack the 
Dirichlet Divisor conjecture. 
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